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We argue that the physics of interacting Kelvin Waves (KWs) is highly non-trivial and cannot 
be understood on the basis of pure dimensional reasoning. A consistent theory of KW turbulence 
in superfluids should be based upon explicit knowledge of their interactions. To achieve this, we 
present a detailed calculation and comprehensive analysis of the interaction coefficients for KW 
turbuelence, thereby, resolving previous mistakes stemming from unaccounted contributions. As a 
first application of this analysis, we derive a new Local Nonlinear (partial differential) Equation. 
This equation is much simpler for analysis and numerical simulations of KWs than the Biot-Savart 
equation, and in contrast to the completely integrable Local Induction Approximation (in which 
the energy exchange between KWs is absent), describes the nonlinear dynamics of KWs. Secondly, 
we show that the previously suggested Kozik-Svistunov energy spectrum for KWs, which has often 
been used in the analysis of experimental and numerical data in superfluid turbulence, is irrelevant, 
because it is based upon an erroneous assumption of the locality of the energy transfer through 
scales. Moreover, we demonstrate the weak non-locality of the inverse cascade spectrum with a 
constant particle-number flux and find resulting logarithmic corrections to this spectrum. 



Q 

U 
d 



(N 
> 

m 
m 



ON 

o 



Physical background, methodology and overview of 
results 



A. Kelvin waves (KWs) in superfluid turbulence 

The role of Kelvin Waves (KWs) in the dissipation of 
energy in zero temperature quantum turbulence has long 
been discussed within the quantum turbulence commu- 
nity. It is widely believed that KWs extend the transfer 
of a constant energy flux from the fully 3D Kolmogorov- 
like turbulence at large scales, through a crossover mech- 
anism at scales comparable to the inter-vortex distance, 
to smaller scales via a local KW cascade on quantized 
vortices. Much theoretical work has been done recently, 
including the conjecture of a power-law scaling for the 
KW cascade made by Kozik and Svistunov in 2004, the 
KS-spectrum 

Nevertheless there remain important unanswered ques- 
tions in quantum turbulence: 

- What are the relative roles of KWs and the other pro- 
cesses, e.g. vortex reconnections, in the transfer of energy 
to small scales? 

- What are the dominant physical mechanisms in the 
classical-quantum crossover range? Two alternative sce- 
narios were put forward for this range: firstly, one relying 
on the idea that the polarization of vortex tangles sup- 
press vortex reconnections, which lead to a bottleneck 
hump 0, [Bj], and the second, implies that reconnections 
play an active role in removing the bottleneck Q. 

- If the KWs do play a key role at small scales, what 
kind of interaction processes are important for the trans- 
fer of energy towards smaller scales? Is it the resonant 
wave- wave interactions, or a linear process of wave num- 
ber evolution due to a large-scale curvature and/or slow 
time dependence of the underlying vortex line, or any 
other possibility? 



- If the KW energy transfer is dominated by the six- wave 
scattering, can one safely assume, as in that this pro- 
cess is local, in the sense that the fc-waves (with a given 
wave vector k) are mainly affected by fc'-waves (with a 
given wave vector k' , where k! = \k'\ is of the same order 
as k = |fe|, with contributions of fe'-waves with k' -C k 
and k' 3> k being vanishingly small? 



In this paper we do not address all of these problems, 
particularly the ones about the role of reconnections and 
about the structure of the crossover range. We restrict 
our attention to the nonlinear interactions of weakly non- 
linear KWs propagating on a single straight vortex line. 
This corresponds to the small-scale range of superfluid 
turbulence at near-zero temperature, where, because of 
the short wavelengths of KWs, one can ignore the influ- 
ence of the neighboring vortex lines within the tangle, 
and assume that nonlincarity, being weak, is still strong 
enough for the nonlinear evolution to proceed faster than 
the large-scale (space and time) changes in the underly- 
ing vortex line. 



Within this idealized setup, our immediate goal is to 
revise and advance the theory of weak-wave turbulence 
of KWs. In particular, we have clarified the structure 
of the nonlinear KW interactions, corrected the theory 
by including previosuly unaccounted leading-order con- 
tributions to the effective wave Hamiltonian, and explic- 
itly calculated the interaction coefficients. Furthermore, 
we have used these results to: 

- Firstly, derive a simple local nonlinear equation (J5J for 
describing KW turbulence. 

- Secondly, to check locality of the KS-spectrum assumed 
in the previous theory. 



2 



B. Hierarchy of the equations of motion for KWs 



H = H 



is defined as follows US 



A commonly accepted model of superfluid turbulence 
comprises a randomly moving tangle of quantized vortex 
lines which can be characterized by the mean intervortex 
distance I and the vortex core radius a <C I. There 
are two approaches in dealing with the vortex core. The 
first one is a "microscopic" model, in which the core is 
resolved: it is based on the Gross-Pitaevski equation, 



at 



0, 



(1) 



where iff is the so-called condensate wave function. 

This model was systematically derived for the Bose- 
Einstein condensate in super-cold atoms, and not for liq- 
uid Helium. Nevertheless, it is frequently used for de- 
scribing superfluid flows in Helium because it contains 
several essential features of such superfluids, i.e. vortex 
quantization, acoustic waves (phonons) in the presence 
of a condensate, and the description of a gradual (non- 
singular) vortex line reconnection. 

However, the Gross-Pitaevski equation can be costly to 
study, and one often resorts to using the so-called Biot- 
Savart formulation of the Euler equations for ideal classi- 
cal fluids, exploiting the fact that far away from the vor- 
tex cores the Gross-Pitaevski dynamics is isomorphic to 
the ideal classical flow via the Madelung transformation. 
In the Biot-Savart model, the vortices are postulated by 
a cutoff in the equations for the vortex line elements. 
Namely, the equations used are @, H, EH : 



k f ds x (r — s) 
4vr / |r - s| 3 : 



(2) 



with a cutoff at the core radius ao, i.e. integration is over 
the range |r — s| > ao. Here k = 2ith/m is the quantum 
of velocity circulation, m is the particle mass. In what 
follows, we will adopt the Biot-Savart equation (BSE) as 
a staring point for our derivation. 

To consider the KW system, one has to start with 
an equilibrium state corresponding to an infinitely long 
straight vortex line and perturb it with small angle dis- 
turbances. This corresponds to a setup of weakly nonlin- 
ear KWs which are dispersive, and that can be decribed 
by weak- wave turbulence theory. For this, one has to 
parametrize the transverse displacement vector (x, y) of 
the perturbed line by the distance z along the unper- 
turbed line (the latter lies along the Cartesian z-axis): 



w(z, t) = x(z, t) + i y(z, t). 



(3) 



In terms of w, the BSE can be written in a Hamiltonian 
form 



. dw _ 5H{w,w*} 
dt Sw* 



(4) 



where d/Sw* is the functional derivative, the asterisk 
stands for the complex conjugation and the Hamiltonian 



BSE K 

H = — 

4-7T 



[l + Re{w'* (z x )w' (z 2 ))] dz 1 dz 2 
\J (zi - z 2 ) 2 + |w(zi) - w{z 2 )\ 2 



(5) 



where we have used the notation w'(z) = dw/dz. 

Then, one must expand in two small parameters: the 
perturbation inclination w' < 1 and Ag -1 <C 1, where 
Ao = \a.{i/ao). Such a simultaneous expansion is not 
easy. This is because the leading order in 1/Ao gives an 
intcgrable model called the Local Induction Approxima- 
tion (LIA) [13, EH with a Hamiltonian 



H 1A = f Ao 



dzy/l + \w'(z)\ 



(6) 



Because of the integrability, LIA conserves an infinite 
number of integrals of motion, and wave resonances are 
absent in all orders, which prevents energy exchange be- 
tween KWs. Thus, LIA appears to be too simple, and to 
describe the KW energy transfer, one has to go to next 
order in 1/Ao. 

The second difficulty is that the lowest order process, 
the four- wave resonances, are absent for such one dimen- 
sional systems with concaved dispersion relations. Thus, 
one must also go to the next order in small w' . The com- 
bination of these two facts makes finding the effective 
interaction Hamiltonian for KWs a hard task. 

For numerical analysis purposes, the BSE model © 
is also quite challenging because of the nonlocal (in the 
physical space) integral that has to be computed, espe- 
cially when one has to resolve a wide range of turbulent 
scales and when the waves are weak, so that the evolution 
times are long. 

Thus, there is a clear need for a simpler model for non- 
linear KWs, which would be local in the physical space 
like LIA (i.e. represent a nonlinear partial differential 
equation), but unlike LIA would be capable of descript- 
ing the energy transfer over turbulent scales. Motivated 



by this need, an ad hoc model was introduced in 13 1 
which has the simplest possible form with all the scal- 
ing properties and solutions of the original BSE model 
preserved. This model was called the Truncated-LIA (or 
TLIA); it has the following Hamiltonian 



TLIA K 

H = 4^° 



1 



w 



dz . 



(7) 



The name TLIA arises from the fact that it can formally 
be obtained by expanding the LIA Hamiltonian in \w'\ 2 
and truncating at the fourth order. This truncation leads 
to the breaking of the LIA integrability, while preserving 
all the important scalings. 

The TLIA model turned out to be very efficient and 
useful for numerical simulations [l3j , even though it was 
suggested ad hoc , and motivated by the need for greater 
simplicity, rather than derived from first principles. In 
the present paper, we will obtain a very simple, local in 
z, nonlinear partial differential equation which is isomor- 
phous to TLIA for weakly nonlinear KWs, and which is 



obtained rigorously from BSE by asymptotic expansions 
in 1/Aq and w' , and by a subsequent identification of the 
dominant type of interaction wave sextets. 



C. The new Local-Nonlinear Equation for Kelvin 

waves 

Our main goal in the paper is, being based on the BSE, 
is to systematically derive an effective motion equation 
for KWs, "as simple as possible, but not more". From 
one side, this new equation will include dynamical inter- 
action (with energy exchange) of KWs (which are absent 
in LIA), and from the other, will contain the leading con- 
tribution obtained from local elements (in contrast to the 
TLIA model). Moreover, at the same time, being drasti- 
cally simplier than the BSE itself. The resulting equation 
reads: 



dw 
~dt 



k d 


[H 


dw 




dw 


Air dz 


dz 


') 


~dz~ 



= 0. 



(8) 



Local Nonlinear Equation (LNE) for KWs. 



The variable w(z,t) is related to w(z,t) via a weakly 
nonlinear canonical transformation of type 



w(z, t) = w(z, t) 
The dimensionless parameter 



0(w' 3 ) 



(9) 



where A is the characteristic wavelength of the the 
KWs. 

It is instructive to represent LNE ((SJ in a Hamiltonian 
form 



■ dw S rr LNE r *l 



(12a) 



where 



^ LNE K 
All 





A 


dw 


2 1 


dw 


6 " 


/ 


dz 


12 


dz 





dz , (12b) 



LNE Hamiltonian. 



In the weakly nonlinear case, the LNE model is isomor- 
phic to the TLIA model ([7]), in which the variable w is 
related to w via the canonical transformation mentioned 
above and a proper reseating. We repeat that the TLIA 
model was introduced in [13[ in an ad-hoc way, requir- 
ing that it must have the same basic scaling properties 
as the original BSE @, but be simpler than the latter. 
In the present paper, the LNE © is derived from the 
BSE systematically and, therefore, it is an asymptotically 
rigorous equation in the relevant KW turbulence limits. 
In spite of the remarkable simplicity of the LNE ©, 
its asymptotical derivation is rather cumbersome, which 
possibly explains the fact that it has never been obtained 
before. To reach our goals we proceeded as follows. 



A = A -7-^- A o-2; 

where 7 = 0.5772 ... is the Euler constant. The replace- 
ment Ao A is equivalent to a replacement of ao by an 
effective vortex core radius a = ao exp(7 + 3/2) ~ 8a , 
in the equation A = ln(f/a). 

We entitled the result © as the Local Nonlinear Equa- 
tion (LNE) to stress its three main features: 

Locality (in physical space) of the KW interactions: ac- 
cording to Eq. ([8]) , the evolution of the KW ampli- 
tude w(z,t) depends only on the slope dw(z,t)/dz 
of KWs and on its curvature (via d 2 w(z, t) / dz 2 ) at 
the same point z. 

Nonlinearity plays a crucial role in the LNE ((SJ, as it is 
responsible for the energy transfer among the KW 
modes. Notice that LIA, being formally a nonlin- 
ear equation, allows (in the framework of the in- 
verse scattering formulation) a linear formulation, 
in which the energy exchange between KWs is ex- 
plicitly absent. 

The word Equation (not an "Approximation" ) stresses 
the fact that LNE ijHJ) is asymptotically exact in the 
triple limit 

w' <1, A _1 <1, and A<l, 



(10) 

D. Challenges of the KW turbulence. Nonlocality 
of the six-wave process. 

The first part of the paper is devoted to an explicit cal- 
culation of the effective six- wave interaction coefficient in 
the limit of small perturbation angles and small values 
of 1 / A that culminate in the derivation of the LNE . 
As we mentioned above, the respective expansion in two 
small parameters is not easy and rather cumbersome. To 
keep our presentation reasonably short and transparent, 
in Sec. II Bl we describe only the main steps of the calcu- 
lation and move details to the Appendices. In this part, 
we fix a set of rather important technical errors made in 
Refs. [H,i], thereby preparing the mathematical problem 
for further analysis. 

As an important application of the obtained result, we 
dedicate Sec. HI] to an analysis of the cascades of energy 
E and particle-number N caused by the 3 O 3-scattering 
of KWs in the context of the KS-conjecture [l|. The KS 
theory predicts a power-law energy spectrum with con- 
stant E-3.uk, E KS (k) (x k~ 7 / 5 , and with constant A-flux, 
E N (k) oc These spectra can only be valid if 

they are local. We show in Sec. Ill CI that the KS spec- 
trum is strongly non-local, so that the KW dynamics 
are dominated by interactions of fc'-waves with fc' <C k. 
(11) Therefore, the KS spectrum is physically irrelevant, i.e. 
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it cannot be realized in Nature. In Sec. Ill CI we also 
demonstrate that the iV-flux spectrum is weakly nonlocal 
in the sense that the dynamics of the k- waves are equally 
effected by all the fc'-waves with k' < k. In Sec. Ill Dt we 
"fix" locality of this spectrum by a logarithmic correc- 
tion, E N (k) oc fc- 3 [ln(fc^)]- 1 / 5 , Eq. pjl . 

Establishing nonlocality of the six- wave theory is a key 
step to the desired effective description. Indeed, it sug- 
gests that the dominant wave sextets must involve modes 
with widely separated wavelengths. We will show that 
in this case, the effective interaction coefficient looks re- 
markably simple: it is proportional to the product of the 
six wave- numbers of the sextet modes, see pip . This im- 
mediately yields the LNE @. In conclusion, we discuss 
prospects of building a theory based upon the LNE (|8|) , 
which should lead to an alternative result to the invalid 
KS spectrum. 



I. HAMILTONIAN DYNAMICS OF KWS 

A. Introduction to the problem 

To derive an effective KW Hamiltonian leading to the 
LNE ©, we first briefly overview the Hamiltonian de- 
scription of KWs initiated in [l[ and further developed 
in 0, [Hj]. The main goal of Sec. II Bl is to start with 
the so-called "bare" Hamiltonian ([5]) for the Biot-Savart 
description of KWs (J2J) and obtain expressions for the 
frequency Eqs. (|2~Tj) . four- and six-KWs interaction coeffi- 
cients Eqs. ([22]) and ([23)) and their 1/A-expansions, which 
will be used in further analysis. Eqs. (J2J), ([22]) and (j2"3")) 
are starting points for further modification of the KW 
description, given in Sec. II CI in which we explore the 
consequences of the fact that non-trivial four- wave inter- 
actions of KWs are prohibited by the conservation laws 
of energy and momentum: 

u)i + uj 2 = w 3 + w 4 , fei + k 2 = k 3 + fc 4 , (13) 

where ujj = w(fcj) is the frequency of the fc^-wave. Only 
the trivial processes with k\ = fc 3 , k 2 — fc 4 , or k\ = fc 4 , 
k 2 = &3 are allowed. 

It is well known (see, e.g. Ref. [l4|) that in the case 
when nonlinear wave processes of the same kind (e.g. 
1 — >■ 2) are forbidden by conservation laws, the terms 
corresponding to this kind of processes can be eliminated 
from the interaction Hamiltonian by a weakly nonlinear 
canonical transformation. A famous example [14[ of this 
procedure comes from a system of gravity waves on water 
surface in which three-wave resonances u\ = uj 2 + w 3 are 
forbidden. Then, by a canonical transformation to a new 
variable b = a + 0(a 2 ), the old Hamiltonian H(a, a*) is 
transformed to a new one, H(b, b*) , where the three- wave 
(cubic) interaction coefficient V { is eliminated at the 
expanse of appearance of an additional contribution to 
the next order term (i.e. four-wave interaction coefficient 




FIG. 1: Examples of contribution of the triple vertices to 
the four-wave Hamiltonian when three-wave resonances are 
forbidden (left) and contributions of the quartet vertices to 
the six-wave Hamiltonian in resonances when four-wave reso- 
nances are forbidden (right). Intermediate virtual waves are 
shown by dash lines. 

if '*) of the type 

V 3 f(vYY/[u5+u 3 -ui] ■ (14) 

One can consider this contribution as a result of the 
second-order perturbation approach in the three-wave 
processes k\ — > k 3 + k 5 and fc 5 + k 2 — > fc 4 , see Fig. [TJ 
left. The virtual wave k$ oscillates with a forced fre- 
quency uj\ — LU3 which is different from its eigenfrequency 
W5. The inequality u>(ki) — w(fc 3 ) ^ w(|fci — fc 3 |) is a 
consequence of the fact that the three-wave processes 
w(fci) — w(fc 3 ) = uj(\ki — fe 3 |) are forbidden. As the result 
the denominator in Eq. (1141) is non-zero and the pertur- 
bation approach leading to Eq. (IT4|) is applicable when 
the waves' amplitudes are small. 

Strictly speaking our problem is different: as we men- 
tioned above, not all 2 -H- 2 processes (fP3|) are forbidden, 
but only the non-trivial ones that lead to energy exchange 
between KWs. Therefore, the use of a weakly nonlinear 
canonical transformation ^ (as suggested in [l| ) should 
be done with extra caution. The transformation ([9]) is 
supposed to eliminate the fourth order terms from the 
BSE-based interaction Hamiltonian by the price of ap- 
pearance of extra contributions to the "full" six- wave in- 
teraction amplitude W^'^, (|23|) . of the following type 
(see Fig. [TJ right): 

rp4.7 ^5,6 

— — , UJ 7 EE oj(\kx + k 2 - fc 4 |) . (15) 

0J- + CJ 4 — LUl — UJ 2 

Here all wave vectors are taken on the six-wave resonant 
manifold 

CJi + CJ 2 + U>3 = ^4 + ^5 + ^6 j (16) 

ki + k 2 + fe 3 = k 4 + k 5 + k 6 . 

The danger is seen from a particular example when 
k\ — » fc 4 , k 2 — > fc 5 and fe 3 — >■ k 6 , so uj-? — > uj 2 , and the de- 
nominator of Eq. (|15p goes to zero, while the numerator 
remains finite. This means, that the perturbation con- 
tribution (lT5j) diverges and this approach becomes ques- 
tionable. 

However a detailed analysis of all contributions of the 
type (|15|) performed a-posteriori and presented in Sec. II CI 
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demonstrates cancelations of diverging terms with oppo- 
site signs such that the resulting "full" six- wave interac- 
tion coefficient remains finite, and the perturbation ap- 
proach (TH))) appears to be eligible. The reason for this 
cancelation is hidden deep within the symmetry of the 
problem, which will not be discussed here. 

Moreover, finding the "full" Hamiltonian is not enough 
for formulating the effective model. Such a model must 
include all contributions in the same 0(1/ A), namely the 
leading order allowing energy transfer in the fe-space. 
The "full" Hamiltonian still contains un-expanded in 
(1/A) expressions for the KW frequencies. Leaving only 
the leading (LIA) contribution in the KW frequency u>, 
as it was done in [J, [|| , leads to a serious omission of an 
important leading order contribution. Indeed, the first 
sub-leading contribution to u shifts the LIA resonant 
manifold, which upsets the integrability. As a result, the 
LIA part of Wi'2'3 yields a contribution to the effective 
model in the leading order. This (previously overlooked) 
contribution will be found and analyzed in Sec. IIP! 

B. "Bare" Hamiltonian dynamics of KWs 

1. Canonical form of the "bare" KW Hamiltonian 

Let us postulate that the motion of a tangle of quan- 
tized vortex lines can be described by the BSE model ([2j , 
and assume that 

A = ln(£/a ) » 1, (17) 

where ao is the vortex core radius. The BSE can be writ- 
ten in the Hamiltonian form Q with Hamiltonian (|5|). 

Without the cut-off, the integral in H , Eq. ([5]), would 
be logarithmically divergent with the dominant contribu- 
tion given by the leading order expansion of the integrand 
in small z\ — z 2 , that corresponds to H , Eq. 

As we have already mentioned, LIA represents a com- 
pletely integrable system and it can be reduced to the 
one-dimensional nonlinear Schrodinger (NLS) equation 
by the Hasimoto transformation [12[. However, it is the 
complete integrability of LIA that makes it insufficient 
for describing the energy cascade and which makes it 
necessary to consider the next order corrections within 
the BSE model. 

For small-amplitude KWs when w'(z) <C 1, we can 
expand the Hamiltonian in powers of w' , see Ap- 
pendix IA II 

H = H 2 + H 4 + H 6 + ... . (18) 

Here we omitted the constant term Hq which does not 
contribute to the equation of motion Eq. (@}. Assuming 
that the boundary conditions are periodical on the length 
C (the limit kC ^ 1 to be taken later) we can use the 
Fourier representation 

w(z, t) = «T 1/2 ^2 a ( fe > exp(ifcz) . (19) 

k 



The bold face notation of the one-dimensional wave vec- 
tor is used for convenience only. Indeed, such a vector 
is just a real number, k € M, in our case for KWs. For 
further convenience, we reserve the normal face notation 
for the length of the one-dimensional wave vector, i.e. 
k = \k\ E In Fourier space, the Hamiltonian equa- 
tion also takes a canonical form: 

. da(k,t) 5H{a,a*} 

1 r = -= — t, — r" ■ (20a) 

dt Sa*(k,t) v ' 

The new Hamiltonian H is the density of the old one: 

H{a, a*} = H{w, w*}/£ = H 2 + H 4 + H 6 + . . . (20b) 

The Hamiltonian 

^ 2 = X! UJk ak a *k ' ( 20c ) 

k 

describes the free propagation of linear KWs with the dis- 
persion law Uk = ui(k), given by Eq. (|21[) . and the canon- 
ical amplitude ak = a(k,t). The interaction Hamiltoni- 
ans H4 and He describe the four- wave processes of 2 O 2 
scattering and the six-wave processes of 3 O 3 scattering 
respectively. Using a short- hand notation a,j = a(kj,t), 
they can be written as follows: 

H4, = - T\'^ aia 2 ala\ , (20d) 

1+2=3+4 

n& = 3^ W if,3 aia 2 a 3 ala* 5 a* 6 . (20e) 

1+2+3=4+5+6 

Here T?;| = T(k 1 ,k 2 \k 3 ,k i ) and Wtf 6 3 = 
W(ki, k 2 , A; 3 1 fc 4 , fc 5 , fe@) are "bare" four- and six-wave 
interaction coefficients, respectively. Summations over 
k\ . . . &4 in H4 and over fei . . . feg in He are constrained 
by ki + k 2 = fe 3 + fc 4 and by fei + k 2 + k 3 = fe 4 + fc 5 + fe 6 , 
respectively. 

2. A-expansion of the bare Hamiltonian function 

As will be seen below, the leading terms in the Hamil- 
tonian functions ojk, T\' 2 and are proportional to 
A, which correspond to the LIA ©. They will be de- 
noted further by a front superscript " A ", e.g. A W/t, etc. 
Because of the complete integrability, there are no dy- 
namics in the LIA. Therefore, the most important terms 
for us will be the ones of zeroth order in A, i.e. the ones 
proportional to A = 0(1). These will be denoted by a 
front superscript " 1 " , e.g. 1 ojk, etc. 

Explicit calculations of the Hamiltonian coefficients 
must be done very carefully, because even a minor mis- 
take in the numerical prefactor can destroy various can- 
celations of large terms in the Hamiltonian coefficients. 
This could change the order of magnitude of the answers 
and the character of their dependence on the wave vectors 
in the asymptotical regimes. Details of these calculations 
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are presented in Appendix IA 11 whereas the results are 
given below. 

Together with Eq. (fTU)) . the Kelvin wave frequency is: 



where 



K A fc 2 / 47T , 

- K fc 2 ln(fc£)/47r 



The "bare" 4-wave interaction coefficient is: 



Tl : 4 2 = A T^ + 1 T^ + 0(A- 1 ), 
A T^2 = -Afcifc 2 fc 3 fc4/47r, 
1 lf$ = - (5fc 1 fc 2 fc 3 fc 4 +^ r i; 2 4 )/f67r 



The function is symmetric with respect to k\ <-} k 2l 
fc,3 <-> k4 and {fci,fc2} <-> {^3,^4}; its definition is given 
in Appendix I A 21 

The "bare" 6-wave interaction coefficient is 



(21a) 
(21b) 
(21c) 



(22a) 
(22b) 

(22c) 



^tli + 'wiii + OiA- 1 ), (23a) 



9A 

87TK 

9 

327TK 



kik 2 k 3 k4k 5 k e 



'7 k^kzkik^ks 



(23b) 



^i 4 ; 2 5 ; 3 6 ) -(23c) 



H{a,a*} 

n 2 



yv 1,2,3 



' 1,2,3 



"1.2 — 



H{b,b*} = H 2 +H A + H 6 
J2 ^ k b k b* k , H± = 0, 



He, — — 



E 

1+2+3=4+5+6 



wtii hb 2 b 3 bib* 5 b* 6 , 



1 



8 ^ 

ij^m— 1 p,q,r— 4 

rpj-.m rpq.p 
r,j-\-m — r i.p-\-q—i 
q r, j+m— r 
j, m 

Wl + w 2 - W 3 - W 4 
A 



(24a) 
(24b) 

(24c) 
(24d) 
(24e) 



rpq,r rpl,J 

m, q-\-r—m p, i+j- 



(24f) 



fi^ + ^^ + ^A- 1 ) . 



The Q-terms in the full six-wave interaction coefficient 



4,5,6 
1,2,3 



can be understood as contributions of two four- 



wave scatterings into resulting six- wave process via a vir- 
tual KW with k = kj + k m — k r in the first term in Q 
and via a KW with k = k q -\-k r — k m in the second term; 
see Fig. [TJ right. 



The function 6^ 1 '2 '3 i s symmetric with respect to ki <H» 
k 2 O fc 3 , ^4 O fes O &6 and {fei, fc 2 , k^} +> {fc4, &5, fee}; 
its definition is given in Appendix IA 31 



Note that the full expressions for Uk, T\ A and W^'I'f 
do not contain £ but rather ln(l/ao). This is natural be- 
cause in the respective expansions £ was introduced as 
an auxiliary parameter facilitating the calculations, and 
it does not necessarily have to coincide with the inter- 
vortex distance. More precisely, we should have used 
some effective intermediate length-scale, £ e g, such that 
£ <g; £ cS <c 2n/k. However, since £ g is artificial and 
would have to drop from the full expressions anyway, we 
chose to simply write £ omitting subscript "eff" . Cance- 
lation of £ is a useful check for verifying the derivations. 



2. 1/A-expansion of the full interaction coefficient W 1 



4,5,6 
2,3 



Similarly to Eq. (|23a[) . we can present W{ 2Z in the 
1 / A-expanded form: 

Wlii = A WlH + 'Win + OiA- 1 ) , (25) 

Due to the complete integrability of the KW system in 
the LIA, even the six-wave dynamics must be absent in 
the interaction coefficient "^W^'j'g. This is true if func- 



tion W 



Aw4,5,6 



1,2,3 



2,3- 

vanishes on the LIA resonant manifold: 



Aw4,5,6 ,4,5,6 
vv 1,2,3 1,2,3 



Ao4,5,6\ _ 



2,3 , 



(26a) 



C. Full "six-KW" Hamiltonian dynamics 



Full six-wave interaction Hamiltonian He 



where 



1-4,5.6 

1,2,3 



A?>4,5,6 
il 1,2,3 



S(ki + k 2 + k 3 - ki - fc 5 - fe 6 ) , (26b) 



— \k 2 
Air [kl 



kr: kc 1 



(26c) 



Importantly, the four-wave dynamics in one- 
dimensional media with concaved dispersion laws 
oj{k) are absent because the conservation laws of en- 
ergy and momentum allow only trivial processes with 
fei = &3J fe2 = fe4, or fei = fe 4 , k 2 — k%. This means 
that by a proper nonlinear canonical transformation 
{a, a*} — > {b, b*}, H4 can be eliminated from the Hamil- 
tonian description. This comes at a price of appearance 
of additional terms in the full interaction Hamiltonian 



Explicit calculation of A W 123 m Appendix [Bl shows that 
this is indeed the case: the contributions A W^'2'3 and 
A Q in Eq. (|24d[) cancel each other (see section 1 in 
Appendix |B|) . (Such cancelation of complicated expres- 
sions was one of the tests of consistency and correctness of 
our calculations and our Mathematica code). The same 
is true for all the higher interaction coefficients: they 
must be zero within LIA. 
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Thus we need to study the first-order correction to the 
LIA, which for the interaction coefficient can be schemat- 
ically represented as follows: 



L W 1 



•4,5,6 



2.:! 



h 1 o 4 ' 5,6 

r l 1 * 1,2,3 



2*9 1 



4,5,6 
2,3 



AO 

x 



' 3 V 1,2,3 ' 
AO ' 



(27a) 
(27b) 
(27c) 



Here 1 W is the A°-order contribution in the bare vertex 
W, given by Eq. (|23c|) . 1 Q is the A°-order contribution 
in Q which consists of \Q and \Q originating from the 
part l T in the four- wave interaction coefficient T, and 
\Q originating from the Kl corrections to the frequencies 
in Eqs. (JHeJ and ((242). Explicit Eqs. jB3| for }<2^3> 

2Q1 2 3 an d 3Q12 3 are presented in Appendix lB 21 They 
are very lengthy and were analyzed using Mathematica, 
see Sec. ITTT21 



is calculated in Appendix IB 31 

W 3,4,5 _ lw3,4,5 , 153,4,5 
yV k,l,2 — VV fc,l,2 + D fe,l,2 ' 



lo4,5,6 _ 2?r 1^4,5,6 
Q l,2,3 — g K "1,2,3 



(^+^)X'2 5 .3 

[kj - h) A 



i={l,2,3} 
j={4,5,6} 



(29a) 
(29b) 



where dj(-) = d(-)/dkj. 

(j2"5|) represents a correct effective model and, will serve 
as a basis for our future analysis of KW dynamics and 
kinetics. However, to make this equation useful we need 
to complete the calculation of the effective interaction co- 
efficient W^' 4, 2 and simplify it to a reasonably tractable 
form. The key for achieving this is in a remarkably simple 
asymptotical behavior of W^' 4 '^, which will be demon- 
strated in the next section. Such asymptotical expres- 
sions for will allow us to establish nonlocality of 
the KS theory, and thereby establish precisely that these 
asymptotical ranges with widely separated scales are the 
most dynamically active, which would lead us to the re- 
markably simple effective model expressed by the LNE 



D. Effective six-KW dynamics 

1 . Effective equation of motion 

The LIA cancelation (|26ap on the full manifold (fT6|) is 
not exact: 

A WIHSIH f(njftt) , 1 ) . 

o 3,4,5 _ , . 

The residual contribution due to 1 O fe ' 1 ' 2 has to be ac- 
counted for - an important fact overlooked in the previ- 
ous KW literature, including the formulation of the effec- 
tive KW dynamics recently presented by KS in 2] . Now 
we are prepared to take another crucial step on the way 
to the effective KW model by replacing the frequency 
ujk by its leading order (LIA) part (I21b|) and simultane- 
ously compensating for the respective shift in the reso- 
nant manifold by correcting the effective vertex T-L. This 
corresponds to the following Hamiltonian equation, 

= A Wfc6fe (28) 

at 

fc+l+2=3+4+5 

where the constraint fc + k\ + fc 2 = k% + + k§ holds. 
Here W^' 4 ' 2 is a corrected interaction coefficient, which 



2. Analysis of the effective interaction coefficient W^'j'j 

Now we will examine the asymptotical properties of 
the interaction coefficient which will be important for 
our study of locality of the KW spectra and formula- 
tion of the LNE (J5J . The effective six-KW interaction 
coefficient W^' 4 ' 2 consists of five contributions given by 
Eqs. (H7J) and @. The explicit form of W^f J involves 
about 2 x I0 4 terms. However its asymptotic expansion in 
various regimes, analyzed by Mathematica demonstrates 
very clear and physical transparent behavior, which we 
will study upon the LIA resonance manifold 

k + fci + fc 2 = fe 3 + fc 4 + fc 5 , (30a) 
fc 2 + k\ + k\ = k\ +k\ + ki . (30b) 

If the smallest wavevector (say fcs) is much smaller than 
the largest wave vector (say fc) we have a remarkably 
simple expression: 

W* 1 a -7^-fcfcifc 2 fc 3 fc4fc5 , (31) 

' ' 47TK 

aa min{fc,fci,fc 2 ,fc 3 ,fc 4 ,fc 5 } 
max{fc, k\ 1 k-2, fc 3 , k±, k§} 

We emphasize that in the expression (|3~Tj) . it is enough 
for the minimal wave vector to be much less than the 
maximum wave number, and not all of the remaining 
five wave numbers in the sextet. This was established 
using Mathematica and Taylor expanding W^' 4 ' 2 with 
respect to one, two and four wave numbers [20|]. All of 
these expansions give the same leading term as in (|31l) . 
see Apps. IB 41 and IB 51 
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The form of expression (|3ip demonstrates a very sim- 
ple physical fact: long KWs (with small fc-vectors) can 
contribute to the energy of a vortex line only when they 
produce curvature. The curvature, in turn, is propor- 
tional to wave amplitude b k and, at a fixed amplitude, 
is inversely proportional to their wave-length, i.e. oc k. 
Therefore, in the effective motion equation each bj has 
to be accompanied by kj, if kj <C k. This statement is 
exactly reflected by formula (|3Tj) . 

Furthermore, a numerical evaluation of W^'^ on a 
set of 2 10 randomly chosen wave numbers, different from 
each other at most by a factor of two, indicate that in 
the majority of cases its values are close to the asymp- 
totical expression (f5T|) (within 40%). Therefore, for most 
purposes we can approximate the effective six-KW inter- 
action coefficient by the simple expression (f3Tj) . Finally, 
our analysis of locality seen later in this paper, indicates 
that the most important wave sextets are those which 
include modes with widely separating wavelengths, i.e. 
precisely those described by the asymptotic formula (I3T1) . 

This leads us to the conclusion that the effective model 
for KW turbulence should use the interaction coefficient 
(|3"TT) . Returning back to the physical space, we thereby 
obtain the desired Local Nonlinear Equation (LNE) for 
KWs given by ©. 

As we mentioned in the beginning of the present paper, 
LNA is very close (isomorphous for small amplitudes) to 
the TLIA modelJT]) introduced and simulated in [13j. It 
was argued in |13| that the TLIA model is a good al- 
ternative to the original Biot-Savart formulation due to 
it dramatically greater simplicity. In the present paper 
we have found a further support for this model, which 
is strengthened by the fact that now it follows from a 
detailed asymptotical analysis, rather than being intro- 
duced ad hoc. 



Partial contributions to the 6-wave effective interaction 
coefficient 



One sees that Eqs. (132]) for the partial contributions 
involve the artificial separation scale £, which cancels out 
from 1 W = X W + \Q + \Q + \Q. This is not surprising 
because the initial expressions Eqs. (f23| do not contain I 
but rather ln(f/ao). This cancelation serves as one more 
independent check of consistency of the entire procedure. 

Notice that in the KS paper contributions (|32dj) 
and (|32e|) were mistakenly not accounted for. Therefore 
the resulting KS expression for the six-wave effective in- 
teraction coefficient depends on the artificial separation 
scale I. This fact was missed in their numerical simula- 
tions [l| . In their recent paper Q , the lack of contribu- 
tion ( |32d[) in the previous work was acknowledged (also 
in [I3|), but the contribution (|32e|) was still missing. 



II. 



KINETIC DESCRIPTION OF KW 
TURBULENCE 



A. Effective Kinetic Equation for KWs 

The statistical description of weakly interacting waves 
can be reached [l4| in terms of the kinetic equation (KE) 
shown below for the continuous limit kC 2> 1, 



dn(k,t)/dt = St(k,t), 



(33a) 



for the spectra n(k,t) which are the simultaneous pair 
correlation functions, defined by 



(b(k,t)b*(k',t)) 



2tt 

T 



6(k-k')n(k,t), 



(33b) 



where (. . .} stands for proper (ensemble, etc.) averaging. 
In the classical limit [2l| . when the occupation numbers 
of Bose particles N(k,t) > I, n(k,t) = HN(k,t), the 
collision integral St(fc, t) can be found in various ways [3, 
SLED], including the Golden Rule of quantum mechanics. 
For the 3 O 3 process of KW scattering, described by 
the motion Eq. 



It would be instructive to demonstrate the relative im- 
portance of different partial contributions, 1 W , \Q, \Q, 
\Q and 1 S [see Eqs. fl27J) and (HH)] to the full effective 
six-wave interaction coefficient. For this, we consider 
the simplest case, when four wave vectors are small, say 
ki,k2,k 3 ,k 5 — > 0. We have (see Appendix IB 5ft : 



W ~ 


■ -1 






(32a) 


\Q 
w 


1 

• + 2 


-\Hki) 


~6 ln ¥' 


(32b) 


\Q 
w 


I 

' + 2 




-IhA 

6 k ' 


(32c) 


IQ 
w 


• +1 


^-Hki) 


1, *3 

+ 6 ln T' 


(32d) 


l s 
w 






6 ln T- 


(32e) 



st 3 « 3 (fc) 



7T 

12 



W 



3,4,5 
k,l,2 



x[n, 



n — no — n A — n 5 



x nA,nin2n37i47i5 dk\ dk2 dk% dk^ dk§ 



(33c) 



KE (|3"3"|) conserves the total number of (quasi)-particles 
N and the total (bare) energy of the system A E, defined 
respectively as follows: 



N 



rife dk . 



co k n k dk 



KE 



has a Rayleigh- Jeans solution, 
T 

n T (k) — —7 ■ , 
h A co k + /i 



(34) 



(35) 



which corresponds to thermodynamic equilibrium of 
KWs with temperature T and chemical potential /i. 
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In various wave systems, including the KWs described 
by KE (|33a[) . there also exist flux-equilibrium solutions, 
n E (k) and n N (fc), with constant fc-space fluxes of energy 
and particles respectively. The corresponding solution for 
n E (fc) was suggested in the KS-paper [l| under an (unver- 
ified) assumption of locality of the 75-flux. In Sec. Ill C[ 
we will analyze this assumption in the framework of the 
derived KE (p£?)) . and will prove that it is wrong. The 
TV- flux solution 7i N (fc) was discussed in 0. in sec, urn 
we will show that this spectrum is marginally nonlocal, 
which means that it can be "fixed" by a logarithmic cor- 
rection. 



B. Phenomenology of the E- and N-flux 
equilibrium solutions for KW turbulence 

Conservation laws (1341 for E and N allow one to intro- 



duce the continuity equations for n k and A E k = A u> k n k 
and their corresponding fluxes in the fc-space, ix k and e k : 



dn k 



at 

d A E k 

at 



+ 



— = , fi k = - J St 3o3 (fc) dk 

de k 



dk 



0, e k 



W St 3o3 (fc) dk 



(36a) 
(36b) 



In scale-invariant systems, when the frequency and in- 
teraction coefficients are homogeneous functions of wave 
vectors, Eqs. (|36]l allow one to guess the scale-invariant 
flux equilibrium solutions of KE (j3"3")l [l4j |: 



n E (fc) = A E k~ 



(37) 



Here A E and A N are some dimensional constants. Scaling 
exponents x N and x E can be found in the case of local- 
ity of the N- and Effluxes, i.e. when the integrals over 
fci, . . . k 5 in Eqs. (l3"6l and (|33c[) converge. In this case, 
the leading contribution to these integrals originate from 
regions where fci ~ fc 2 ~ fc 3 ~ fc 4 ~ k$ ~ fc and thus, the 
fluxes (|36p can be estimated as follows: 



fc 5 [W(fc,fc,fc|fc,fc,fc)] 2 n^(fc). 



(38a) 
(38b) 



Stationarity of solutions of Eqs. (1361) require constancy 
of their respective fluxes: i.e. fi k and e k should be fc- 
independent. Together with Eqs. (|3"5|) this allows one to 
find the scaling exponents in Eq. (|37p . 

Our formulation (|33l) of KW kinetics belongs to the 
scale- invariant class [22| : 

A tjfe oc fc 2 , and for V 7/ 
W^fc, 7/fci , r?fe 2 | ?7 fc 3, 77^4, ^5) 
= T 1 6 W(k,k 1 ,k 2 \k 3 ,k i ,k 5 ) . 

Estimating W(fc, fc, k\k, fc, fc) ~ fc 6 /«; and A oj k — ftAfc 2 in 
Eqs. (p8| . one gets for iV-flux spectrum [3j: 



n N (k) ~ (/x/t/A) fc 3 , x N = 3 , 



(39a) 



and for 75- flux KS-spectrum [l|: 

n E (fc) ~ (eK 2 ) 1/5 fc" 17/5 , x E = 17/5 . (39b) 



Non-locality of the N- and iS-fluxes by 
3 -o- 3-scattering 



B for KWs with 



Consider the 3 O 3 collision term 
the interaction amplitude Wj^ 33. Note that in (|33c[) 
J (ifej are one-dimensional integrals rffcj ■ Let us ex- 
amine the "infrared" (IR) region ( fcs <C fc, fci, fc2, fc 3 , fc4 ) 
in the integral (I33c[) , taking into account the asymp- 
totics (l3~Tj) , and observing that the expression 



r 3,4,5 r/Ao 3,4, 5\ / -1 , -1 , 



n 9 — 71, — n 1 



x n k n\ 7i 2 n 3 774 77,5 

x nfc ni n 2 n 3 n 4 n 5 — n 5 ~ fc 5 _;7: . 

Thus the integral over k$ in the IR region can be factor- 
ized and written as follows: 



;^fc 2 7i 



(fc 5 ) dk 



■5 oc 



i/t 



k 2 - 



C dh 



(40) 



The factor 2 here originates from the symmetry of the 
integration area and evenness of the integrand: f_ = 
2 J °°. The lower limit in this expression should be 
replaced by the smallest wave number where the assumed 
scaling behavior (|37[) holds, and moreover, it depends 
on the particular way the wave system is forced. For 
example, this cutoff wave number could be l/£, where £ is 
the mean inter-vortex separation £, at which one expects 
a cutoff of the wave spectrum. The crucial assumption 
of locality, under which both the E-tiux (KS) and the 
TV-flux spectra were obtained, implies that the integral 
(j40| is independent of this cutoff in the limit £ 0. 
Clearly, integral (j40|) depends on the IR-cutoff if x > 3, 
which is the case for both the 75-flux (KS) and the TV-flux 
spectra (f39|). Note that all other integrals over fci, fc 2 , fc 3 
and &4 in p3cj) diverge exactly in the same manner as 
the integral over fcs, i.e. each of them leads to expression 

Even stronger IR divergence occurs when two wave 
numbers on the same side of the sextet (e.g. fci and fc 2 , 
or fc 3 and k^, etc.) are small. In this case, integrations 
over both of the small wave numbers will lead to the 
same contribution, namely integral (|40j) . i.e. the result 
will be the integral (j40]l squared. This appears to be the 
strongest IR singularity, and the resulting behavior of the 
collision integral Eq. (133c[) is 




2-x 



dkn 



(41) 



When two wave numbers from the opposite sides of the 
sextet (e.g. k 2 and ks) tend to zero simultaneously, we 
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get an extra small factor in the integrand because in this 
case (n^ 1 + n^ 1 — n^ 1 — x ) — > 0. As a result we 
get IR convergence in this range. One can also show IR 
convergence when two wave numbers from one side and 
one on the other side of the sextet are small (the resulting 
integral is IR convergent for x < 9/2). 

Divergence of integrals in Eq. (|33cp means that both 
spectra (|39p with a; N = 3 and x E = 17/5 > 3, obtained 
under opposite assumption of the convergence of these 
integrals in the limit £ — > oc are not solutions of the 
3 <H> 3-KE (|33c[) and thus cannot be realized in Nature. 
One should find another self-consistent solution of this 
KE. Note, that the proof of divergence at the IR lim- 
its is sufficient for discarding the spectra under the test, 
whereas proving convergence would require considering 
all the singular limits including the ultra-violet (UV) 
ranges. However, we have examined these limits, too. 
At the UV end we have obtained convergence for the KS 
and for the inverse cascade spectra. Thus the most dan- 
gerous singularity appears to be in the IR range, when 
two wave numbers from the same side of the wave sextet 
are small simultaneously. 



D. Logarithmic corrections for the iV-fiux 
spectrum (|39al) 

Note that for the A^-flux spectrum (|39a[) n N (fc) oc fc~ 3 , 
that the integrals (|4"0")) and (|4T|) diverge only logarithmi- 
cally. The same situation happens, e.g. for the direct en- 
strophy cascade in two-dimensional turbulence: dimen- 
sional reasoning leads to the Kraichnan-1967 [15| turbu- 
lent energy spectrum 



E(k) oc fc 



-3 



(42a) 



for which the integral for the enstrophy flux diverges log- 
arithmically. Using a simple argument of constancy of 
the enstrophy flux, Kraichnan suggested [l6[ a logarith- 
mic correction to the spectra 



E{k) oc k- 3 \n~ 1/3 (kl), 



(42b) 



that permits the enstrophy flux to be fc-independent. 
Here I is the enstrophy pumping scale. 

Using the same arguments, we can substitute in 
Eq. (|33cl) . a logarithmically corrected spectrum n N (k) oc 
k~ 3 \n~ x (k£) and find x by the requirement that the re- 
sulting iV-flux, /ifc, Eq. (I36a[) will be fc-independent. Hav- 
ing in mind that according to Eq. (I38a[) ^ oc n^, we can 
guess that x — 1/5. Then, the divergent integral (j40|) will 
be oc ln 4 / 5 (fc£), while the remaining convergent integrals 
in Eq. (|33cl) will be oc ln~ 4 / 5 (fc£). Therefore, the resulting 
flux [ik will be fc-independent as it should be [Hj]. So, our 
prediction is that instead of a non-local spectrum (|39a[) 
we have a slightly steeper log-corrected spectrum 



n N (fc) 



.1/5 



fc 3 ln 1/5 (fc^) 



(43) 



The difference is not large, but the underlying physics 
must be correct; as one says on the Odessa market: "We 
can argue the price, but the weight must be correct" . 



Conclusions 

In this paper, we have derived an effective theory of 
KW turbulence based on asymptotic expansions of the 
Biot-Savart model in powers of small 1/A and small non- 
linearity, by applying a canonical transformation elim- 
inating non-resonant low-order (quadric) interactions, 
and by using the standard Wave Turbulence approach 
based on random phases [l3 |. In doing so, we have fixed 
errors arising from the previous derivations, particularly 
the latest one by KS Q, by taking into account pre- 
viously omitted and important contributions to the ef- 
fective six- wave interaction coefficient. We have exam- 
ined the resulting six-wave interaction coefficient in sev- 
eral asymptotic limits when one or several wave numbers 
are in the IR range. These limits are summarized in a 
remarkably simple expression pip . This allowed us to 
achieve three goals: 

• To derive a simple effective model for KW turbu- 
lence expressed in the Local Nonlinear Equation 
([HJ . In addition to small 1/A and the weak nonlin- 
earity, this model relies on the fact that our findings 
show, for dynamically relevant wave sextets, the in- 
teraction coefficient is a simple product of the six 
wave numbers, Eq. (|3ip . For weak nonlinearities, 
the LNE is isomorphic to the previously suggested 
TLIA model [lj]. 

• To examine the locality of the i?-flux (KS) and the 
A^-flux spectra. We found that the KS spectrum is 
non-local and therefore cannot be realized in Na- 
ture. 

• The A^-flux spectrum is found to be marginally non- 
local and could be "rescued" by a logarithmic cor- 
rection, which we constructed following a qualita- 
tive Kraichnan approach. However, it remains to 
be seen if such a spectrum can be realized in Quan- 
tum Turbulence, because, as it was shown in [l7j . 
the vortex line reconnections can generate only the 
forward cascade and not the inverse one (i.e. the re- 
connections produce an effectively large-scale wave 
forcing) . 

Finally we will discuss the numerical studies of KW 
turbulence. The earliest numerics by KS were reported in 
[lj. They claimed that they observed the KS spectrum. 
At the same time they gave a value of the E-Rux con- 
stant ~ 10~ 5 which is unusually small. We have already 
mentioned that this work failed to take into account sev- 
eral important contributions to the effective interaction 
coefficient, and thus these numerical results cannot be 
trusted. In particular, we showed that their interaction 
coefficient must have contained a spurious dependence on 
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the scale I which makes the numerical results arbitrary 
and dependent on the choice of such a cutoff. In addi- 
tion, even if the interaction coefficient was correct, the 
Monte-Carlo method used by KS is a rather dangerous 
tool when one deals with slowly divergent integrals (in 



this case J Q 



-7/5 



dx). 



On the other hand, recent numerical simulations of 
the TLIA model also reported agreement with the KS 
scaling (as well as an agreement with the inverse cascade 
scaling) (l3| . How can one explain this now when we 
showed analytically that the KS spectrum is non-local? 
It turns out that the correct KW spectrum, which takes 
into account the non-local interactions with long KW's, 
has an index which is close (but not equal) to the KS 
index, and it is also consistent with the data of [l3j]. We 
will report these results in a separate publication. 
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Appendix A: Bare interactions 

1. Actual calculation of the bare interaction 
coefficients 

The geometrical constraint of a small amplitude per- 
turbation can be expressed in terms of a parameter 



This allows one to expand Hamiltonian ([5]) in powers of 
e and to re-write it in terms consisting of the number of 
wave interactions, according to Eq. (TTH1) . KS found the 
exact expressions for H 2 , H4 and H§ [1]: 



H 2 
Hi 
H 6 



k f dz\dz 2 
8n J \zi-z 2 \ 

k f dz\dzi 
327T J \zi-z 2 \ 

k f dz\dz 2 
64W \ Zl -z 2 \ 



2Re 
lie 



(w *{z 1 )w (z 2 ) 
4 -Ae 2 Re (w '* (zi)w \z 2 ) 
6e 4 i?e (w *(zi)w (z 2 )^j — 5c 



,(A2) 



e(zi,z 2 ) = \w(z\) - w(z 2 )\/\z 1 - z 2 \ < 1. 



(Al) 



The explicit calculation of these integrals was ana- 
lytically done in [l3j], by evaluating the terms in (|A2[) 
in Fourier space, and then expressing each integral as 
various cosine expressions [l[. Hamiltonian (|A2I) can 
be expressed in terms of a wave representation variable 
afe = a(k,t) by applying a Fourier transform (|19p in the 
variables z\ and z 2 , (for details see [l|, HH). The result 
is given by Eqs. (|20[) . in which the cosine expressions for 
ujk , and W 4 23 were done in [l[ . In our notations they 



arc 



uo k = ±[A-B\, T 34 = ±-[6D-E\, W 4 ^ = jL[3P- 5Q] , where (A3) 

ZTT 47T 47W 

A=f°°—k 2 C k 1 B=r^[l-C k ], D= r^i-d-a-CP-CP + Ci + CU + CZ], 

J an z ~ J 0,0 z — J 00 z — 

r°° dz 

E = — [hki (C 4 + d~ C 43 - C7 2 4 ) + hk 3 (C 3 + Ci - C 43 - Cl) + k 3 k 2 (C 3 + C 2 - C 43 - cl) 

+ k t k 2 (C 4 + C 2 - C 43 - Cf ) ] , (A4) 

r°° dz 

P = -^k 6 k 2 [C 2 - Cf - C 23 + Cl 3 - Cl + Cf + C 4 3 - Cl + C 6 - C 56 - Cl + Cl 6 - C 46 + C 45e + Cl 6 - C12] , 

J an z - 



an 

dz- 



Q = I - c 4 ~ c\ + ct - c 6 + c 46 + cf - ct e ~c 5 + c 45 + c\ - c 45 + c 65 - c 456 - c? 6 + c 23 

Jan z - 

- O3 + O3 + O i3 ~ U 13 + U 3 — U 3 - 13 + U 2 + U 3 — U 3 - 13 + U 2 — U 3 + O12 + 2 — U 2 J 

I 

Here the variable, z~ — \zi — z 2 \ and the expressions C, are cosine functions such that C\ = cos(fciz_), C\ 
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cos((fc 4 - kt)z-), Cf = cos((fe 4 + k 5 - fci)z_), Cff = 
cos((fc 4 + k 5 — ki — k 2 )z^) and so on. The lower limit 
of integration ao is the induced cutoff of the vortex core 
radius ao < \z\ — z%\. 

The trick used for explicit calculation of the analyt- 
ical form of these integrals was suggested and used in 
[l3j |. First one should integrate by parts all the co- 
sine integrals, so they can be expressed in the form of 
cos(z) 

dz. Then, one can use a cosine identity for 



this integral [1J 



cos(z) , , , , f a ° cos(z) -1 

—dz = -7 - ln(oo) - / — dz (A5) 



00 (_ a 2\ k 
. 7 _ ln(ao) _g__ 



- 7 -ln(|a |) + O(ag), 



where 7 = 0.5772 ... is the Euler Constant. Therefore, 
in the limit of a small vortex core radius ao, we can ne- 
glect terms of order ~ a^ and higher. For example, let's 
consider the following general cosine expression that can 

z -3 cos(/Cz)dz, where JC is 

an expression that involves a linear combination of wave 
numbers, i.e. K. = k\ — fc 4 . Therefore, integration by 
parts will yield the following result for this integral: 



be found in Eqs. (IA4|) : 



cos(JCz) 



dz = 



ao 

K 

1 



cos(/Cz) 


OO 


K sin(/Cz) 


2z 2 


+ 

a 


2z 



cos(JCz) 



dz 



cos(/Cao) /Csin(/Ca ) 



2al 



2a n 



2 



cos(y) 



fCao 



v 



dy 



We then expand cos(/Cao) and sin(/Cao) in powers of ao, 
and apply the cosine formula (|A5|) for the last integral, 
where in the last step we have also changed integration 
variables: y — K,z. The final expression is then 



cos(/Cz) , 1 3/C 2 K 2 , , ... m . 2 . 



Explicit equations for T\ 2 and 5 4 ff are given below in 
Appendices IA 21 and IA 31 In the main text we introduced 
A = A - 7 - 3/2. Writing A = \n(£/a), we see that 
a = a e 7+3 / 2 ~ 8a . 



2. Bare 4-wave interaction function T 1 , 2 

A rather cumbersome calculation, presented above, re- 
sults in an explicit equation for the 4-wave interaction 
function J 7 12 in Eqs. (|A6|) and (|22bl) . Function J-f'2 is a 

,3,4 



symmetrical version of F J'^: ^12 = l^i'^j where the 
operator {. . . } g stands for the symmetrization k\ k 2l 
k% <r> ki and {ki,k 2 } f-> {fc3,fc 4 }. In its turn F 3 ; 4 . is 
defined as following: 



F 3 ; 4 = ]T /C 4 ln(|/CK) 
/ce/Ci 

+2^ ^ fc^j /C 2 ln(|/C|^) 



(A7a) 



The j denotes sum of four terms with = 
{(4, 1), (3, 1), (3,2), (4,2)}, K, is either a single wave vec- 
tor or linear combination of wave- vectors that belong to 
one of the following sets: 





= {- 


[1] 


rb] 


r[ 3 ]r[V [t], + [ 4 V g]}> 


£41 


- V 


[ 4 ] 


+ [1] 


,-n-m , (A?b) 


/C31 


= V 


[ 3 ] 




r [ 43 ]r [1]} , 


IC32 


= V 


[ 3 ] 


+ w 


r [ 43 ]r [?]} , 




- V 


[ 4 ] 




r[ 43 ]r[ 4 ]} . 



Here we used the following shorthand notations with 
a,/3,7 = 1,2,3,4: [«] = fc Q , [9] = -fe^ , g] = 



/3 



y ] = + ^7 j [/3-y] — 



/3 



fc 7 , and 



or ~~ signs before [. . . ] should be understood as prefac- 
tors + 1 or — 1 in the corresponding term in the sum. For 
example: 



By applying a similar procedure to the other cosine in- 
tegrals, we find that all terms of negative powers of ao, 
(that will diverge in the limit a — > 0) actually cancel in 
the final expression for each interaction coefficient. Ap- 
plying this strategy to all interaction cofficients, we get 
the following analytical evaluation of the Hamiltonian 
functions 131: 



/C 4 ln(|/C|<?) for/Ce{~[i]}is -fc 4 ln(M), 
/C 4 ln(|/C|£) for fC € { + H}} is + (fc 4 - k 2 ] 
kikj K? \ai\K\l) for i = 4, 3 = 1, 
/CG{"[ 43 ]}is -fc 4 fci(fe4 + fc3) 4 ln(|fe 4 4 



4 ln(|fc 4 -fc 2 |£), 
■k 3 \£) . 



1 12 



W- 



456 
123 



Kk 2 r 

47T 
1 

loTT 
9 

327TK 



A = ln(£/a ) , 
Ao - 7 - I - Hki) 



(A6) 



fcife 2 fc 3 fc4(l + 4 7 - 4A ) - 

fcifc 2 fc 3 fc 4 fc 5 fc6(i - 4 7 + 4A ) - gt;^ 



3. Bare 6- wave interaction function Q j_' 2 * s 



Function 5 4 '2'3 = ^ ^ 123 j" • The operator {. . . } g 
stands for the symmetrization ki ■n- k 2 O fc3, k± ■n- 
fcs ■(-> fee a nd {fci,fc2,fc3} ^ {^4,^5,^617 and G 4 '^ is 
defined as following: 
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-o4,5,6 
1,2,3 



fc 6 fc 2 /C 4 ln(|/C|£) 

K.3 

l - J2 >C 6 H\m , (A8a) 



K£K 2 
1 

+ 18 



where = w i + ^2 — ^3 — u>4- We want to 

compute this equation on the LIA manifold (f3"U)) . To do 
this we express two wave vectors in terms of the other 
four [23[ using the LIA manifold constraint (f3"U)) : 



where 



Ka = 



/C4 = 



|2ji L23J 

r56i — r6i +[56i 
I 1 ' [3J> L3 J ^ 



{-I 



+ r5 



r4i +r45i + r4 



r46i 



r456i 



r6i +r6i 



L23h L2h L2 J j L23h LlJ> 



[46 1 
[3 J- 



[12]}, (A8b) 



+ r45i +r5i -r45i +r65i -r456i -r56i +r 1 -r 
[ J; llJi [l J ; I J> [ J) [l J: [23J, I 

"""[t]' + [l3], [lsl^IjiL [3 6 ] > [13]' + 



},-[t% + [i2), + [t\,-[2}}. (A8c) 



fcl 



(fc 3 - fe) (fc 2 - fc 3 ) 

fc + &2 — &3 — &5 
(fc 3 - fc) (fc 2 - fc 3 ) 

fe + k 2 - k 3 - fc 5 



fc 5 , (B2a) 
k + k 2 -k 3 . (B2b) 



Appendix B: Effective six-KW interaction coefficient 

1. Absence of 6- wave dynamics in LIA 

According to Eqs. (|24d[) and (|24e[) . the expression for 
A ^i2'3 * s gi ven by 



A n 4,5,6 
V 1,2,3 



A p,q,r 



1 3 

- y 



^ K 1,2,3 ^ 1,2,3 ' 
3 6 

EA p,<J,r 
2,j,m— 1 p 7 q,r= 4 

Arpj,m Arpq,p 

r, j+m— r i,p+g— i 

j, m 

At-t<?, r Arpi, j 

m, g+r — m p, i+j—p 



AO™ 1 ' ' 



m, q-\-r— m 



(Bla) 
(Bib) 



(Blc) 



Then A W^x'a'3 ^ s easnv simplified to zero with the help of 
Mathematica. This gives an independent verification of 
the validity of our initial Eqs. (l24l) for full interaction cof- 



ficient A W / i' 2 ' 3 which is needed for the calculations of the 



O(l) contribution 1 W 123 - Another way to see the can- 
celation is to use the Zakharov-Schulman variables [l8| 
that parameterise the LIA manifold (|30l) . 

2. Exact expression for 1 W 



We get expressions for \Q, \Q and 3 Q, introduced by 
Eqs. (|2T|) . from Eqs. (|24djl and pie]) . Namely: 



l n 4,5,6 
1^1,2,3 



3 6 

\ E E 

i,j,m— 1 p.q,r— 4 



r, j-\-m- 



lrpq,p 

i, p+q- 



A^j+m- 



Arpq,r lrpl,J 

m, q+r—m p, i+j- 

AO m ' 9+ r— m 



(B3a) 



i4,5,6 
'l,2,3 



3 6 

;E E 

i,j,m— 1 p,g,r= 4 



1 jii, 



r, j-\-k—r 



Arpq,p 

i, p+q- 



AurJ+m- 
3,m 



Irpq, r 

tci. q+r- 



Arpl, J 

P, i+j— P 



An m ' 9+' r - 



(B3b) 



£,4,5,6 
^1,2,3 



3 

: E 

i^j^m— 1 3 



6 

E 

,q,r= 4 



AtX), m 

r, j-\-m- 



Arpq,p 

i,p+q- 



l n r, j+m— r 
J, ™ 



Arpq,r 

m, q+r- 



( Ar> m -< 

y ^ L q. r 



Arpi,j 

P, i+3-P . l^m,q+r- 



m, q+r — m ' 



(B3c) 



Again, using Mathematica we substitute Eqs. (|B2[) into too cumbersome to be presented here. But we will ana- 
Eqs. (|B3a|) - (|B3cl) . Clearly, the resulting equations are 
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lyze them in various limiting cases, see below. 



This is the required expression Eq. (|29b[) . 



3. Derivation of Eq. {29b]) for 



First of all, let us find a parametrization for the full 
resonant manifold, by calculation of the correction to the 
LIA parametrization (|B2[) . namely 



fci = A ki + % 



fcj 



A fc 4 



where A fci and A fc4 are given by the right-hand sides of 
Eqs. (|B2|) respectively. Corrections : fci and x k 4 are found 
so that the resonances in fc, Eq. (|30ap . and (full) u> are 
satisfied. The resonances in fc fixes x fei = ^4. Then the 
w-resonance in the leading order in 1 /A gives 



4,5,6 
1,2,3 



L k 



d\j 4 
dfc 4 



Thus 



^4 



l-K "1,2,3 



An (fe 4 — fci) 



(B5) 



(B6) 



This allows us to write down the contribution of A W from 
the deviation of the LIA resonant surface: 



1 £4,5,6 _ 1, 
J l,2,3 — K l 



Aw4,5,6 
2,3 



d A W? 



dki 

2tt ^4.5.6 <A 

An ^ 



- ^4 „, 1,2,3 

afci 

+ gl) A ^l 4 2 5 3 
(fc 4 - fci) 



©(A" 1 ) 

(B7) 



with dj(-) = dj(-)/dkj. It is obvious that instead of fci 
and fc 4 we could use parametrizations in terms of other 
pairs ki and kj with i = 1, 2 or 3 and j = 4, 5 or 6. This 
enables us to write a fully symmetric expression for 1 S: 



1 £4,5,6 _ 

'-'1,2,3 — 



2VT ^ 4 ,5,6 



9Ak 



2.3 



E 



2,3 



i={l,2,3} 
j={4,5,6} 



(kj ki) 



(B8) 



4. Analytical expression for W on the LIA 
manifold when two wave numbers are small 



Let us put together the coefficients to the interaction 
(B4) coefficient Wjj^j given in {27a]) . (f27b|) . ([27c]) . ([23c]) and 



I, and use in these expressions the formulae obtained 
in the previous appendices and the parametrization of 
the LIA surface (|B2j) . Using Mathematica, and Taylor 
expanding W^'^ with respect to one wave number, e.g. 
fcs, we obtain a remarkably simple result, - expression 

Now we will consider the asymptotical limit when two 
of the wave numbers, say k 2 and k§ (let them be on 
the opposite sides of the resonance conditions) , are much 
less than the other wave numbers in the sextet. Using 
Mathematica and Taylor expanding W^'^ with respect 
to two wave numbers fc2 and k$, we have 



lim W 



3,4,5 

fc,l,2 



-- — k 2 k 2 kjk 5 

4-7TK 



(B9) 



Simultaneously, see Eq. (|B2I) : 



lim fci — > k 3 , 
fc 2 ->■ 
fc 5 ->• 



lim fc4 — » fc 

fc 2 -> 
fc 5 -> 



(BIO) 



Therefore, (|B9|) coincides with (I3T1) . Note that this was 
not obvious a priori, because formally pip was obtained 
when fcs is much less than the rest of the wave numbers, 
including fc2. 

For reference, we provide expressions for the different 
contributions to the interaction coefficient Wf'^'o given 
in Eqs. {27]) and {29]). For fc 2 , fc 5 -> 0: 
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3Q -> 



3 

47TK 



47TK 



47TK 



3 

47TK 



k 2 k 2 k 2 k 5 



k 2 k 2 k^k 5 



k 2 k 2 k 2 k 5 



k 2 k 2 k 2 k 5 



--^—k 2 k 2 k 2 k 5 
Attk 



ln(M)-^(49 



2x(12 + 7ar) In |a; 
3 



(1 -x) 2 (7+ 10a; + 7a: 2 



7^ mil 



•In II 



ln(fc£) 



2 

2 (9a; 2 
3 



1 

48 



59 (l-^) 2 (9 + 10^ + 9x 2 ) ln 



m(fci) 



14a; - 6 - 
X -7. 



In lie j 



(1 + .x) 



1 — x' ' ' a; 
(1 - a) 2 (l + a: 2 ) 



mil 



'0 



In |1 



2 1 - 5x + a:3 mbi 



(1 + xf 



1 1 



1-a 

x 



In |1 



6 1 - x 



1 1 
1- - 



In 



a 



6 1 — x 



In la 



x = fea/fe 



(Blla) 
(Bllb) 

(Bile) 
(Blld) 

(Bile) 



r 



Another possibility is for two small wave numbers to be 
on the same side of the sextet. We have checked that on 
the resonant manifold, this also leads to p"T]) . 



5. Analytical expression for W on the LI A 
manifold when four wave numbers are small 

Now let us, using Mathematica, calculate the asymp- 
totic behavior of W when four wave vectors are smaller 
than the other two; on the LIA manifold this automat- 
ically simplifies to ki,k 2 ,k 3 ,k^ <C k, fej (remember that 
on the LIA manifold k\ and k^ are expressed in terms of 
the other wave numbers using Eq. (|B2|) . thus from (IBlOp ) 
we have 



lim W 



3,4,5 

fc,l,2 



Attk 



k 2 k 2 k 2 k 5 . (B12) 



Again, we have got an expression which coincides with 
(1311) . We emphasize that this was not obvious a priori, 
because formally pip was obtained when fcs is much less 
than the rest of the wave numbers, including k\,k 2 ,k 3 . 

Therefore we conclude that the expression pip is valid 
when &5 is much less than just one other wave number in 



the sextet, say k, and not only when it is much less than 
all of the remaining wave numbers. 

For a reference, we give the term by term results for 
the limit k\,k 2 ,k 3 ,k^ <C k,k±: 



l w 


->• 


3 

-- — k 2 k 2 klk 5 
Attk 


-1 - 




+ , 


\Q 


-> 


3 

-- — k 2 k 2 k\k b 
Ann 


1 

+ 2 


--ln(^) 


1, V 

-6 ln ¥ 


\Q 


->• 


3 

-- — k 2 k 2 k\k b 
Ann 


1 

+ 2 




1, k 3 ~ 


IQ 


-> 


3 

-- — k 2 k 2 klk 5 
Attk 


+1- 




f 6 ln T. 


x s 


-> 


3 

- k k 2 k<]k§ 
Attk 





+ 


1, k 3 - 



The sum of this contributions is very simple: 



x W-> k 2 k 2 k 2 k 5 \+l 

Attk j l 
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